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A new explanation for the longstanding puzzle of the tiny branching fraction of 
?// —> (tk relative to that for J/ip —> pn was proposed. In the case of J/tp, we argue 
that this decay is dominated by a higher Fock state in which the cc pair is in a 
color-octet 3 Si state and annihilates via the process cc —> qq. In the case of the 
ijj', we argue that the probability for the cc pair in this higher Fock state to be 
close enough to annihilate is suppressed by a dynamical effect related to the small 
energy gap between the mass of the ip' and the DD threshold. 


1 Introduction 

Usually, a well-defined puzzle of physics provides us with a gateway to look at 
something new. We get insight into the physics behind the puzzle through seek¬ 
ing its solution. There is such a long-standing puzzle in charmonium physics 
which is called as the “p — tt puzzle” of J/^> and ip' decays. The puzzle is 
that, J/ip decays to pir but ip 1 does not, which contradicts our understanding 
of the decay processes. The discrepancy between the conventional theoreti¬ 
cal expectation and the experiments is greater than 65. This extremely large 
number clearly tells us that there is something unusual in charmonium decays. 
We hope to learn something from this 15 year mystery. Recently, Chen and 
Braaten proposed a new explanation for it El. That’s what I am going to talk 
about. 

J/if) and ip 1 are nonrelativistic bound states of a charm quark and its 
antiquark. Their decays into light hadrons are believed to be dominated by 
the annihilation of the cc pair into three gluons. In order to annihilate, the c 
and c must have a separation of order l/m c , which is much smaller than the 
size of the charmonium state. Thus the annihilation amplitude for an S-wave 
state like J/ip or ip 1 must be proportional to the wavefunction at the origin, 
ip( r = 0). The width for decay into any specific final state h consisting of light 
hadrons is therefore proportional to |'0(O)| 2 . The width for decay into e + e _ is 
also proportional to |^(0)| 2 . This leads to the simple prediction that the ratio 
of the branching fractions for ip' and J/ip is given by the “15% rule”: 

_ Q* a = (>«±2.3)%. a) 
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This rule follows simply from the fact that ip' can make a transition to J/ip 
and Xc by emitting 7 r’s and 7 . This decreases the branching ratio of the ip' —> h 
annihilation decay. This rule is obeyed by most of the J /ip and ip' annihilation 
decay modes. The p- 7r puzzle is that the prediction (1) is severely violated in 
the pn and several other decaymhannels. The Mark II collaboration found that 
Q P k < 0.6% and Qk*k < 2%l3. Recent data from the BES collaboration has 
made the puzzle even sharpen. They obtained Q pn < 0.23 %, Qk*+k- < 0.64%, 
and Q k ,oko = (1.7 ± O.6)%0. Thus the suppression of Q relative to the 15% 
rule is about 10 for K*°K° and greater than 65 for pir. 


2 Previous explanations 

The p — 7 r puzzle challenges our conventional theory. It has attracted broad 
interest. In literature, several different explanations have been njjoposed. A 
summary of those explanations has recently been given by Chaoa. Hou and 
SoniB suggested that J/ip — * pn is enhanced by a mixing .of the J/ip with a 
glueball O that decays to pn. Brodsky, Lepage, and TuanB emphasized that 
J/ip —> pn violates the helicity selection rule of perturbative QCD, and argued 
that the data requires O to be narrow and nearly degenerate with the J/ip. 
Present data from BES constrains the mass and width of thej-glueball to the 
ranges | mo — mj/^,\ < 80 MeV and 4 MeV < To < 50 MeVll This mass is 
about 700 MeV lower than the lightest J PC = 1" glueball observed in lattice 
simulations of QCD without dynamical quarks □. There are explanations of 
the p-n puzzle that involve the dependence of the decay amplitude on the 
energy of the charmonium state. Karl and Robertscl suggested that the decay 
proceeds through cc — > qq followed by the fragmentation of the qq into pn. 
They argued that the fragmentation probability is an oscillatory function of 
the energy-which could have a minimum near the mass of the ip'. Chaichian and 
Tornqvist Ed pointed out that the suppression of ip' decays could be explained 
if the form factors for two-body decays fall exponentially with the energy as 
in the nonrelativistic quark model. There are other explanations of the p-n 
puzzle that rely on the fact that, there is a node in the radial wavefunction 
for ip' , but not for J/ip. Pinskyliil suggested that this node makesjf// —> pn a 
“hindered Ml transition” like J/ip —> ^. 7 . Brodsky and KarlinerEd suggested 
that the decay into pn proceeds through intrinsic charm components of the p 
and n wavefunctions. They argued that the cc pair in the | udcc) Fock state 
of the p + or n + has a nodeless radial wavefunction which gives it a larger 
overlap with J/ip than ip'. Finally, Li, Bugg, and ZouEj have suggested that 
final-state interactions involving the rescattering of a\p and a- 2 p into pn could 
be important and might interfere destructively in the case of the ip'. 
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3 A new explanation 


Recently Braaten and I proposed a new explanation of the p-n puzzle El. We 
argue that the decay J/ip —► pn is dominated by a higher Fock state of the 
J/xl i in which the cc is in a color-octet 3 Si state. The cc pair in this Fock state 
can annihilate via cc —> g* —> qq. The amplitude for forming pn is dominated 
by the endpoint of the meson wavefunctions, where a single q or q carries most 
of the momentum of the meson. The suppression of ip' —> pn is attributed to 
a suppression of the cc wavefunction at the origin for the higher Fock state of 
the tl)'. Such a suppression can arise from a dynamical effect associated with 
the small energy gap between the mass of the il>' and the DD threshold. 


3.1 J/ip and ip' inclusive annihilation decay rate 

It is convenient to analyze charmonium decays using nonrelativistic QCD 
(NRQCD)liil, the effective field theory obtained by integrating out the energy 
scale of the charm quark mass m c . In Coulomb gauge, a charmonium state has 
a Fock state decomposition in which the probability of each Fock state scales 
as a definite power of v, the typical relative velocity of the charm quark. The 
dominant Fock state of the J/ip and the ip' is |cci( 3 Si)), whose probability V 
is of order 1. We denote the color state of the cc pair by a subscript (1 for 
color-singlet, 8 for color-octet) and we put the angular momentum quantum 
numbers in parentheses. The Fock states |cc 8 ( 3 P/) + S), where S represents 
dynamical gluons or light quark-antiquark pairs with energies of order m c v or 
less, have probability V ~ v 2 . The next most important Fock states include 
|cc 8 ( 1 S' 0 ) + S) with V ~ v 3 (or V ~ v 4 if perturbation theory is sufficiently 
accurate at the scale m c v) and |cc 8 ( 3 »Si) + S) with V ~ u 4 . 

The decay of the J/ip into light hadrons proceeds via the annihilation of 
the c and c, which can occur through any of the Fock states. This is expressed 
in the NRQCD factorization formula for the inclusive decay rate li-3, which 
includes the terms 


J/ifj—tl.h. 
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where a s = a s (m c ). The matrix elements are expectation values in the J/ip of 
local gauge-invariant NRQCD operators that measure the inclusive probability 
of finding a cc in the J/ip at the same point and in the color and angular- 
momentum state specified. The matrix element of the cci( 3 <Si) term in (|) 
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is proportional to the square of the wavefunction at the origin and scales as 
v 3 . Its coefficient includes a term of order a 3 from cc —> ggg and a term of 
order a 2 from the electromagnetic annihilation process cc —> 7* —► qq. The 
color-octet terms in ([|) represent contributions from higher Fock states. Their 
matrix elements scale like v 6 , v 7 and v 7 , respectively. Their coefficients are 
all of order a 2 and come from cc —> gg for the ccs( 1 5'o) and ccs( 3 Po) terms 
and from cc —* g* —> qq for the ccs( 3 S'i) term. Note that the coefficients of 
the color-octet matrix elements are two orders of magnitude larger than that 
of (Oi( 3 Si))j/^, which suggests that the higher Fock states may play a more 
important role in annihilation decays than is commonly believed. 


3.2 Color-octet contribution to J/tj} —> pit 

According to a general QCD factorization formula, the amplitude for a two- 
body annihilation decay can be expressed in terms of hard-scattering factors T 
that involve only the scale m c , initial-state factors that involve scales of order 
m c v and lower, and final-state factors that involve only the scale Aqcd- If 
m c was asymptotically large, the dominant terms in the factorization formula 
would have the minimal number of partons involved in the hard scattering. 
Terms involving additional soft partons in the initial state are suppressed by 
powers of v. Terms involving additional hard partons in the final state am sup¬ 
pressed by powers of Aqcd/ m c . As pointed out by Brodsky and Lepagetd, the 
leading term in the asymptotic factorization formula is strongly constrained 
by the vector character of the QCD interaction between quarks and gluons. 
It vanishes unless the sum of the helicities of the mesons is zero. Rotational 
symmetry then requires the angular distribution to be 1 — cos 2 9. This helic- 
ity selection rule is violated by the decay J/if) —> pn. Parity and rotational 
symmetry require the helicity of the p to be ±1 and the angular distribution 
to be 1 + cos 2 9. Since the helicity of the pion is 0, the helicity selection rule 
is violated. Thus the amplitude for J/ip —> pit is suppressed by a factor of 
Aqcd /m c relative to that for generic mesons. The leading contribution is a 
term with a hard scattering factor of order a 2 from the process cci —► uddug. 

Since the charm quark mass m c is less than an order of magnitude larger 
than Aqcd, there can be large corrections to the asymptotic decay amplitude. 
In particular, there can be significant regions of phase space in which some of 
the gluons involved in the hard scattering are relatively soft. It might therefore 
be more appropriate to absorb them into the initial-state or final-state factors. 
In this case, not all of the soft partons in S need be involved in the hard 
scattering, and not all the partons that form the p and 7 r need be produced by 
the hard scattering. For example, there can be a contribution from the Fock 
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state |ccg( 3 S'i) + S) that involves the hard-scattering process cc —► qq, where 
q = u or d. This produces a state | qq + S) that consists of the soft partons 
S together with a q and a q that are back-to-back and whose momenta are 
approximately m c . Such a state has a nonzero overlap with the final state 
\pn). The overlap comes from the endpoints of the meson wavefunctions, in 
which most of the momenta of the p and n are carried by the q and q. When S 
is a qq pair, the p and n mesons are formed from the 2-quark antiquark pairs, 
with the hard quark pair carries most of the momenta of the p and n. 

This contribution to the T-matrix element can be written schematically in 
the form 

^( c cs( 3 5 , i) -*• qq) ^2{p^\qq +S) (ccsCSx) + s\j/i/j). (3) 

99 S 

Note that the initial-state and final-state factors on the right side of (© cannot 
be separated, because they are connected by the sum over soft modes S. The 
T-matrix element (jsj) would contribute to the cc8( 3 Si) term in the factorization 
formula d) for inclusive decays. 

The endpoint contribution in (| leads to a definite angular distribution. 
Since the q and q carry most of the momenta of the mesons, the angular 
distribution of the mesons will follow that of the q and q , which is 1 + cos 2 0. 
Thus (^]) will contribute most strongly to form factors which allow the angular 
distribution 1 + cos 2 9. It will also contribute most strongly to decays into 
mesons like p and n for which most of the momentum can be carried by a single 
q or q. There are also endpoint contributions involving the hard-scattering 
process ccs( 3 Pj) —> gg , which produces a pair of hard gluons with the angular 
distribution 2 + cos 2 9, and ccg( 1 5'o) — 1 > gg, for which the angular distribution 
is isotropic. Their contributions to J/'i/’ —♦ pn are suppressed by the small 
probabilities for most of the momentum of the p or n to be carried by a single 
gluon and by the mismatch between the angular distribution of the gluons and 
the 1 + cos 2 9 distribution of the pn. 

We argue that the color-octet term in (|j) may actually dominate the decay 
rate for J/ip —> pn. We compare the various factors in that term with those 
in the asymptotic amplitude. The hard-scattering factor T in the color-octet 
term in (|||) is only of order a s , compared to a 2 for the asymptotic ampli¬ 
tude. The suppression from the initial-state factor in the color-octet term in 
(H), including the sum over S. might be as little as a factor of v 2 relative 
to the asymptotic amplitude. This follows from the fact that the color-octet 
amplitude contributes in quadrature to the (Os( 3 Si)) term in (|j), which is sup¬ 
pressed by v 4 . As for the final-state factors, (|l|) is suppressed by the endpoints 
of the meson wavefunctions, while the asymptotic amplitude is suppressed by 
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Aqcd / m c from the violation of the helicity selection rule. Considering the 
various suppression factors, it is certainly plausible that the cc8( 3 Si) term in 
(3j) could dominate over the leading contribution from the cci( 3 Si) Fock state. 


3.3 Suppression of ip' —> pn 

We have argued that the decay J/ip —> pn may be dominated by the annihila¬ 
tion of the cc pair in the |cc8( 3 Si) + S) Fock state via cc —> qq. If this is true, 
then the p —w puzzle can be explained by a suppression of this decay mechanism 
in the case of the if'. This suppression can arise from the initial-state factor 
(ccs( 3 iS'i) + S\J/ip) if the cc wavefunction for the |ccs( 3 5'i) + S) Fock state is 
suppressed in the region in which the separation of the cc is less than or of 
order l/m c . Note that it does not require a suppression of the probability for 
the higher Fock state, but just a shift in the probability away from the region 
in which the c and c are close enough to annihilate. A possible mechanism 
for this suppression is a dynamical effect related to the small energy gap be¬ 
tween the mass,of the ip' and the DD threshold. In the Born-Oppenheimer 
approximation^, the cc pair in the dominant |cci) Fock state moves adiabat- 
ically in response to a potential V\ (R) given by the minimal energy of QCD 
in the presence of a color-singlet cc pair with fixed separation R. Similarly, 
the cc pair in the |ccs( 3 5'i) + S) Fock state moves adiabatically in response 
to a potential V$(R) given by the minimal energy of the soft modes S in the 
presence of a color-octet cc pair with fixed separation R. At short distances, 
this potential approaches a repulsive Coulomb potential a s /(6R). At long dis¬ 
tances, the minimal energy state consists of D and D mesons separated by a 
distance R , and V$(R) therefore approaches a constant 2 (Mp — m c ) equal to 
the energy of the DD threshold. A charmonium state spends most of the time 
on the color-singlet adiabatic surface, but it occasionally makes a transition 
to the color-octet adiabatic surface. Since the J/ip is 640 MeV below the DD 
threshold, a cc pair on the color-octet adiabatic surface is far off the energy 
shell. The time spent by the J/ip on this surface is too short for the cc pair to 
respond to the repulsive short-distance potential. Since the wavefunction for 
the |cci) Fock state peaks at the origin, the cc wavefunction for the |ccs + S) 
Fock state should have significant support near the origin. However the mass 
of the ip' is only 43 MeV below D + D~ threshold and 53 MeV below D°D° 
threshold. A cc pair on the color-octet adiabatic surface can be very close to 
the energy shell. The ip' can therefore spend a sufficiently long time on this 
surface for the cc pair to respond to the repulsive short-distance potential. 
This response can lead to a significant suppression of the cc wavefunction at 
the origin for the |ccs + S) Fock state. 
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This suppression effect can be illustrated further as follows. The repulsive 
interactions between c and c in the color-octet higher Fock state |cc 8 + S'} leads 
to the increase of the distance between c and c. This suppresses the annihilation 
probability of the cc pair. Although this suppression effect is complicated we 
can expect that it is sensitive to the lifetime of the color-octet higher Fock state 
in J/ip and ip'. The longer the lifetime is, the larger the suppression effects are. 
In the case of J/■ip, the lifetime of this state is shorter because it is far off-shell. 
An extreme case is the decay of the J/ip or T to 3 hard gluons. After one or 
two hard gluons are emitted, the core quark is far off-shell and the lifetime 
of the intermediate state is very short. The effect of the repulsive interaction 
can be accounted for as part of the radiative correction to the short-distance 
coefficients. Consequently, we need only one parameter ip( 0) to describe the 
annihilation of the cc pair in this process. In the case of ip', the lifetime of 
the intermediate higher Fock state is much longer because it’s very close to 
the DD threshold. The repulsive interactions separate c and c to a very large 
distance until it is screened by a DD state, in which the annihilation rate is 
very tiny. In this case, the amplitude will be enhanced by a factor 1/A E, the 
energy difference between the mass of the ip' and the DD state, which arises 
from the time integral from 0 to oo in perturbative calculation. However, the 
smaller phase space may compensate this enhancement. This is similar to 
the electromagnetic transitions between two different energies in atoms, where 
transition rates are smaller for those with smaller energy difference. 

If the initial-state factor (cc 8 ( 3 <Si) + S| J/ip) in (||) is suppressed for all soft 
partons S , the suppression can be expressed in the form of a relation between 
the NRQCD matrix elements in (j2j): 

(C> 8 ( 3 5i))^< „ (e > 8 ( 3 5 , i)) J/v , 

This inequality can be tested by calculating the matrix elements using Monte 
Carlo simulations of lattice NRQCD. Since the ip' is so close to the DD thresh¬ 
old, it would be essential to include dynamical light quarks in the simulations. 


4 Signatures 

Our proposal implies interesting signatures. It leads to predictions for the 
flavor-dependence of the suppression of the decays of ip' into vector/pseudoscalar 
final states. It can explain the J PC dependence of the suppression of the two- 
body decays of ip'. Also it gives prediction for the angular distribution of these 
two-body decays. 
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/.I Flavor dependence 


One of the signatures of our explanation is the prediction for the flavor- 
dependence of the suppression of the decays o£_7>' into vector/pseudoscalar 
final states. Bramon, Escribano, and Scadron EJ have analyzed the decays 
J/xp VP assuming that the decay amplitude is the sum of a flavor-connected 
amplitude g, a flavor-disconnected amplitude rg , and an isospiii- violating am¬ 
plitude e. The expressions for the amplitudes are given in Ref. lill. They allowed 
for violations of 517(3) flavor symmetry through parameters s and x — 1. The 
authors give two sets of parameters that fit the existing data, one with x = 1 
and the other with x = 0.64. Both sets have e comparable in magnitude to rg 
and an order of magnitude smaller than g. If the decay J/xp —> pn is domi¬ 
nated by endpoint contributions, we can identify g and e with the two terms in 
(j|). While rg may also have endpoint contributions from cc —> gg, we assume 
for simplicity that it is dominated by subasymptotic contributions from the 
cci( 3 5i) Fock state. The amplitudes for the decays xp' —> VP can be expressed 
in a similar way in terms of amplitudes g' , e', and (rg)'. Our explanation of 
the p— 7r puzzle implies that \g'\ is much smaller than \g\, and that e' and (rg)' 
differ from e and rg only by the factor required by the 15% rule. The unknown 
amplitude g' is constrained by the BES data on xp' —► pn and xp' —* K*°K°. 
The upper bound on B(xp' —> pir) gives an upper bound on \g' + e'| 2 , which 
implies that g' lies in a circle in the complex g'-plane. The BES measurement 
of B(xp' —» K*°K°) gives an allowed range for |(1 — s)g' — (1 + x)e'\ 2 , which 
constrains g' to an annulus. The intersection of the interior of the circle with 
the annulus is the allowed region for g'. By varying g' over that- region and 
taking into account the uncertainties in the parameters of Ref. Eli, we obtain 
the predictions for Qvp in Table 1. Measurements of the xp' branching frac¬ 
tions consistent with these predictions would imply that the suppression of the 
vector/pseudoscalar decays is due to the suppression of g'. This would lend 
support to our explanation of the p -7r puzzle. 

There is a possibility that the isospin-violating amplitude e is dominated 
by the contribution from the cc color-singlet higher Fock state |cci( 3 5i) + S) 
of the J/xp and the xp' through an annihilation of the cc —> 7 * —> qq. If this is 
the case, then we expect there to be no suppression of the 15% rule for the xp' 
decay from this amplitude, since it arises from the cc color-singlet state. The 
ratio of the cross sections of e + e _ —» utk and e + e _ —> p + p~ near the J/xp and 
the xp' resonances should then be significantly smaller than the corresponding 
branching ratios at the resonances. This possibility can also be tested by 
observing the decay rates of J/xp , xp' —► e + e~ 7 T 7 r or J/xp , xp' — > p + p~Tnr 
with soft 7r’s, since these processes receive contribution from the |cci( 3 5i) + 5) 


Table I. Predictions for Qvp in units of 1% for all the vector/pseudoscalar final 
states. The values for pn and K*°K° + c.c. were used as input. The columns 
labelled x = 1 and x = 0.64 correspond to the two parameter sets of Ref. 113 . 


VP 

X = 1 

x = 0.64 

P7T 

0 - 0.25 

0 - 0.25 

K*°I<° + c.c. 

1.2 - 3.0 

1.2 - 3.0 

K*+K~ + c.c. 

0 - 0.36 

0 - 0.52 

LOT] 

0 - 1.6 

0 - 1.6 

LOT]' 

10 - 51 

12 - 55 


0.8 - 3.6 

0.4 - 3.0 


0.7- 2.5 

0.5 - 2.2 

PV 

14 - 22 

14 - 22 

pv' 

12 - 20 

13-21 

UJ7T 

11 - 17 

11 - 17 


higher Fock states. 

4-2 J PC dependence 

A solution to the p-n puzzle should also be able to explain the pattern of 
suppression for various J PC states with the same flavors. A preliminary mea¬ 
surement of The axial-vector/pseudoscalar decay mode ip' —► bin by the BES 
collaboration^ gives Qb^ = (24 ± 7)%, consistent with no suppression relative 
to the 15%_rule. A preliminary measurement of the vector/tensor decay mode 
ip' —> pa -2 ell gives Q pa2 = (2.9 ± 1.6)%, which, though suppressed relative to 
the 15% rule, is an order of magnitude larger than Q p7T . This pattern can be 
explained by also taking into account the orbital-angnlar-momentum selection 
rule for exclusive amplitudes in perturbative QCD 1 3. The decay modes &i7r 
and pct 2 both have form factors that are allowed by the helicity selection rule. 
They also both have form factors that violate the helicity selection rule, but 
are compatible with an endpoint contribution from cc —> qq. However, in the 
case of b±n, the endpoint contribution is further suppressed by the violation 
of the orbital-angular-momentum selection rule. Thus we should expect no 
suppression of ip' —■> 6i7t and only a partial suppression of ip' —> pa- 2 - 

We also expect that the decays of the J/ip and the ip' into the isospin 
violating axial-vector/pseudoscalar modes bip and bip' etc. should be seriously 
suppressed relative to &i7r decay modes, because they are dominated by a 
contribution from the leading Fock state which is suppressed by (a e m/ q s ) . 
The decays of the ip' into isospin violating vector/tensor modes such as a\u> 
and fip should not be suppressed relative to the J/ip decays, because they are 
dominated by the contribution from the color-singlet Fock state. A thorough 
analysis of J /ip and ip' decays into axial-vector/pseudoscalar and vector/tensor 
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final states will be presented elsewhere. 


4-3 Angular distribution 

Our proposal also has implications for the angular distributions of two-body 
decay inodes. In general, the angular distribution must have the form 1 + 
a cos 2 8, with — 1 < a < +1. Our solution to the p —n puzzle is based on the 
suppression of a contribution to ip' decays that gives the angular distribution 
1 + cos 2 9. In the case of J/ip —» pir, the angular momentum selection rule 
implies that a = +1. Since J/ip —> bin and ip' —> bin are dominated by the 
contribution from the cc leading Fock state, we expect that a is close to —1. 
The parameter a for the ip' —> pa 2 should be less than a for the corresponding 
J/ip decay. 

5 Conclusion 

In conclusion, we have proposed a new explanation of the p— n puzzle. We 
suggest that the decay J/ip —> pn is dominated by a Fock state in which the 
cc is in a color-octet 3 Si state which decays via cc —> qq. The suppression of 
this decay mode for the ip' is attributed to a dynamical effect that suppresses 
the cc wavefunction at the origin for Fock states that contain a color-octet cc 
pair. Our explanation for the p-n puzzle can be tested by studying the flavor 
dependence of the two-body decay modes of the J/ip and ip', their angular 
distributions, and their dependence on the J PC quantum numbers of the final- 
state mesons. 
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